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ABSTRACT 

Future X-ray observations of galaxy clusters by high spectral resolution missions will 
provide spatially resolved measurements of the energy and width for the brightest 
emission lines in the intracluster medium (ICM) spectrum. In this paper we discuss 
various ways of using these high resolution data to constrain velocity power spectrum 
in galaxy clusters. We argue that variations of these quantities with the projected 
distance R in cool core clusters contain important information on the velocity field 
length scales (i.e. the size of energy-containing eddies) in the ICM. The effective length 
^off along the line of sight, which provides dominant contribution to the line flux, 
increases with R, allowing one to probe the amplitude of the velocity variations at 
different spatial scales. In particular, we show that the width of the line as a function of 
R is closely linked to the structure function of the 3D velocity field. Yet another easily 
obtainable proxy of the velocity field length scales is the ratio of the amplitude of the 
projected velocity field (line energy) variations to the dispersion of the velocity along 
the line of sight (line width). Finally the projected velocity field can be easily converted 
into 3D velocity field, especially for clusters like Coma with an extended flat core in the 
surface brightness. Under assumption of a homogeneous isotropic Gaussian 3D velocity 
field we derived simple expressions relating the power spectrum of the 3D velocity 
field (or structure function) and the observables. We illustrate the sensitivity of these 
proxies to changes in the characteristics of the power spectrum for a simple isothermal 
/3-model of a cluster. The uncertainties in the observables, caused by stochastic nature 
of the velocity field, are estimated by making multiple realizations of the random 
Gaussian velocity field and evaluating the scatter in observables. If large scale motions 
are present in the ICM these uncertainties may dominate the statistical errors of line 
width and shift measurements. 

Key words: X-rays: galaxies: clusters - Galaxies: clusters: intracluster medium - 
Turbulence - Line: profiles - Methods: analytical - Methods: numerical 



1 INTRODUCTION 

Properties of gas motions in the hot intracluster medium 
(ICM) are still little known. It is believed that turbulent 
motions are driven when matter accretes onto the filaments 
or during shocks in the hot gas. Turbulence transfer kinetic 
energy injected on large scales L ~ Mpc to small (unknown) 
dissipative scales I. These two scales are connected with a 
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cascade of kinetic energy, which occurs ove r inertial range 
(|KoImogorovlll94ll : iLandau fc Lifshitj|l966l ). 

Knowing the properties of gas motions in clusters, we 
would be able to address a number of important question, 
e.g., what is the bias in the cluster mass measurements based 
on hydrostatic equilibrium and whether the bias is due to 
the motio ns alone or d ue to the clumping in gas density 
(see, e.g.. iRasia etal.1 120 06: Nagai, Vikhlinin, & Kravt sovl 
I2OO7I : Ijeltema et al.1 I2OO8I : iLau. Kravtsov. fc Nagail 1200*9^ 
what is the ICM turbulent heating rate in clusters (e.g. 
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Figure 1. Simulated spectrum of the Perseus cluster core (cen- 
tral 130 kpc), calculated assuming different Mach numbers of the 
turbulence. At each point along the line of sight the locally emit- 
ted spectrum is calculated using Astrophysical Plasma Emission 
Code (APEC. ISmith et aT] | |200ll 'l') model using known radial den- 
sity and temperature dependencies. The lines are then broadened 
by the thermal and turbulent motions using eq. I I24II (see i|7.2ll 
and the whole spectrum is projected along the line of sight. The 
case of pure thermal broadening of lines is shown in black. This 
plot shows that (i) turbulent broadening is important even for the 
modest level of turbulence and (ii) linos overlap with each other, 
especially if turbulent motions are present. 



IChurazov et alj 12008 ) and what is the r ole of gas mo- 
tions in particle accelera tion (see, e.g. iBrunettil 120061 : 
iBrunetti fc Lazarianlboill ). 

Properties of turbulence in galaxy cluster s were stud- 
ied by means of numerica l simu l ations (e.g., Dolag et al.l 



19991 : 



2005'; 'C assano fc Brunc tti" '20051: iNorman fc BrvanI 

lapichino et al.l l201ll : IVazza ct al.l 1201 ll ). However despite 
the good "global" agreement between all simulations, 
the results on turbulent motions are still controversial, 
mainly due to insufficient resolution of simulations and, 
in particular, low Reynolds number (effective Re < 
1000 in cosmolo gical simulations) and other numerical is- 



sues (see, e.g., Kitsionas et al.l l2009l : lOobler et al.l l2003l : 



iBeresnvak fc Lazarianll2009l ') 



Current generation of X-ray observatories cannot pro- 
vide robust direct measurements of turbulence in the 
ICM. Only XMM RGS grating can provide weak up- 
per limits on velocity amplitude in cool core clus- 
ters (jSanders. Fabian, fc Smithll201ll ). Indirect indications 
of the ICM turbulence come fr om measurements of 
the resonant scatt ering effect (e.g. IChurazov et al.l |2004| : 
I Werner et al.l 120091 ') . from measurements of pressure fluctu- 



ations ( Schucc ker et al. 200J) or surface brightness fluctua- 
tions (C hurazov et al.ll2012f l 

Future X-ray observatories, such as Asiro-iu and 
ATHENj^ with their high energy resolution will allow us to 
measure shifts and broadening of individual lines in spectra 
of galaxy clusters with high accuracy. Combination of di- 
rect measurements of velocity amplitudes with indirect mea- 
surements via resonan t scattering will give us constraints on 
anisotropy of motions (jZhuravleva et al.ll201lh . X-ray polari- 
metric measurements can also provide in formation on gas 
motio ns perpendicular to the line of sight (|Zhuravleva et al.l 
l2010t ). 

Here we discuss the possibility of getting information 
about the length scales of gas motions (e.g. the size of 
energy-containing eddies). We discuss various ways to con- 
strain structure function and power spectrum of gas motions 
via measurement of the projected velocity (shift of the line 
centroid) and the velocity dispersion (broadening of the line) 
as a function of projected distance from the cluster center. 
These ideas are illustrated with a very simple model of a 
galaxy cluster. An application of our methods to simulated 
galaxy clusters will be considered in future work. 

A similar problem of obtaining the structure func- 
tion of turbulence from spectral observations has been 
addressed in the studies of Galactic interstellar medium. 
In particular, it was shown that the width of molecu- 
lar spectral lines increases with the size of a cloud (see 
e.g. [Mvers et al.1 1 19781 : iHever fc BruntI l2004l : iHever et all 
|2009|). This correlation w as interpreted in terms of tur- 
bulent velocity spectrum (JLarsonI Il98ll ). A way to con- 
strain structure function of turbulence in the Interstel- 
lar Medium (ISM) by means of the velocity centroids 
(project ed mean velo c ity) m easur ements was fl rst consid- 
ered by Ivon Hoernerl (Il95ll) a nd iMiinchI () 19581 ) (see also 
iKleiner fc Dickmarj|l983l . ll985l ). Currently several different 
flavors of the velocity centroids method are used for studies 
of the ISM turbulence (see e.g. lEsauivel et al.ll2007h . More 
advanced techniques, such as Velocity Channel Analysis 
(VCA) and Velocity Coordi n ate S p ectrum (VCS) (see e.g. 
Lazarian fc Poeosvanl |2000|, 120081 : IChepurnov fc LazarianI 



20091 ). were developed and applied to t he ISM data in the 



Milky Way and othe r galaxies (see e.g. iPad oan ot al 



IStanimirovic fc Lazarian . 200ll: IChepurnov et al. 201(|). Few 
other methods were also used, am ong them are the Spec- 
tral Correlation Function (SCF) (JRosolowskv et all Il999l : 
IPadoan. Goodman, fc Juvelall2003n and th e Principal Com- 
ponent Analvsis (PCA) (|Brunt et al.ll2003h . 

The ISM turbulence is often s upersonic and compress- 
ible (e.g. lElmegreen fc Scald IJOoi ) . This leads to (i) large 
shifts in the energy centroid of the line compared to the 
thermal broadening and (ii) large amplitude of the gas den- 
sity fluctuations. At the same time individual lines, which 
are used to study the ISM turbulence (e.g. 21 cm line of 
HI or CO lines) are often well separated from other emis- 
sion lines. The regions under study often have very irregular 
structure on a variety of spatial scales. The analysis there- 
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Figure 2. Maps of observables. Left: surface brightness in line (arbitrary units); middle : emissivity-weighted projected velocity field 
in km/s (centroid shift of the line); right: velocity dispersion in km/s (line broadening). Here we consider only the velocity component 
along the line of sight. We assume /3-model of the emissivity with /3 = 0.6 and re = 10 kpc and a random, uniform and homogeneous 
Gaussian velocity field with cored power 3D power spectrum (eq. [CTJ with km = 0.005 kpc^^ and Kolmogorov slope a = —11/3 (see 
Section 2). The size of each map is IMpcxlMpc. 



fore is usually concentrated on separation of the velocity and 
density fluctuations in the observed data, while the thermal 
broadening can often be neglected. 

In contrast, in galaxy clusters the gas motions are 
mostly subsonic. The detection of the gas motions is still 
possible, since we deal with the emission lines of ions of 
heavy elements like e.g. Fe, Ca or S. The atomic weights of 
these elements are large (e.g. for Fe it is 56) and this drives 
pure thermal broadening of lines down (see Fig. [l] and Sec- 
tion l7.2p . The brightest lines in spectra of galaxy clusters are 
often very close to each other. For example, in the vicinity of 
the He-like iron line at 6.7 keV there are forbidden and in- 
tercombination lines and a number of satellite lines, energy 
separation between which is of the order few tens eV (Fig. 
[l]). Density of clusters often have a regular radial structure 
with relatively small amplitude of stochastic density fluctu- 
ations. Analysis of X-ray surface brightness fluctuations in 
Coma cl uster shows that dens ity fluctuations are ~ 2 — 10 
per cent (jChurazov et al.ll2012l ) . Also hydrodynamical simu- 
lations of cluster formati on predict very small clumping fac- 
tors within rsoo (se e e.g. iMathiesen. Evrard. fc Mohiill999l : 
iNagai fc Laull201ll ). Therefore, to the first approximation, 
the contribution of density fluctuation in galaxy clusters can 
be neglected (see Section [7.41 for details), while global ra- 
dial dependence has to be taken into account (Section |4]). 
Another characteristic feature of X-ray observations is the 
importance of the Poisson noise, related to the counting 
statistics of X-ray photons. If one deals with the clusters 
outskirts the high energy resolution spectra will be domi- 
nated by the Poisson noise even for large area future tele- 
scopes, like e.g. ATHENAq. Finally, one can mention, that 
the effects of self-absorptions can potentially be important 
in clusters. Galaxy clusters are transparent in X-rays in con- 
tinuum and in most of the lines. However, some strong lines 
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can have an optical depth ~ few units. Therefore if one mea- 
sures the width of the optically thick line distortions due to 
res onant scattering effect should be taken into account (see 
e.g. IChurazov et al.ll2010l: IWerner et al.ll2009l : Ide Plaa et all 
l2012l : IZhuravleva et al.ll201ll ). 

Presence of several closely spaced emission lines, modest 
level of turbulence (i.e. modest ratio of the turbulent and 
thermal broadenings) , lack of very strong stochastic density 
fluctuations on top of a regular radial structure, and often 
strong level of Poisson noise affect the choice of the simplest 
viable approach to relate future observables and most basic 
characteristics of the ICM gas velocity field. 

The fact that lines in spectra of galaxy clusters are 
very close to each other (e.g. Fig. [T} and the line ratio 
is temperature dependent can be circumvented by esti- 
mating the mean shift and broadening with direct fitting 
of the projected spectra with the plasma emission model 
(possible multi-temperature model). While the small ther- 
mal broadening of lines from heavy elements helps to ex- 
tend the ap plicability of VCA/VCS techniques into subsonic 
regime (see Esguivel et al.ll2003l : iLazarian fc Esquivelll2003l : 
IChepurnov et al.M2010l ). limited spectral resolution of the 
next generation of X-ray bolometers (e.g. FWHM ~ 5eV for 
ASTRO-H) reduces the measured inertial range in the veloc- 
ity domain. Direct application of SCF and PCA methods to 
galaxy clusters can also be challenging, especially when the 
spectra are dominated by the Poisson noise. These problems 
should be alleviated with missions like ATHENA, having 
very large effective area and an energy resolution of ~ 2.5 
eV. 

Below we suggest to use the simplest "centroids and 
broadening" approach as a first step in studying the ICM 
turbulence. This approach assumes that at any given po- 
sition one fits the observed spectrum with a model of an 
optically thin plasma (including all emission lines) and de- 
termines the velocity centroid and the line broadening. This 
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reduces the whole complexity of X-ray spectra down to two 
numbers - shift and broadening. Simple analysis of exist- 
ing hydrodynamic simulations of galaxy clusters shows that 
this approximation does a good (although not perfect) job 
in describing the profiles o f emission lines (see e.g. Fig. 2 in 
llnogamov fc Sunvaevl2003l ) . At the same time this approach 
is the most effective in reducing the Poison noise in the raw 
measured spectra. We show below that in spite of its sim- 
plicity this approach provided an easy way to characterize 
the most basic properties of the ICM velocity field. 

Clearly, more sop histicated metho ds developed for the 
ISM turbulence (e.g. iLazarianl 12003 ') will eventually be 
adopted to the specific characteristics of the ICM turbu- 
lence and X-ray spectra, potentially providing a more com- 
prehensive description of the ICM turbulence, once the data 
of sufficient quality become available. 

The structure of the paper is as follows. In Section 2 
we describe and justify models and assumptions we use in 
our analysis. In Section 3 we specify observables which can 
be potentially measured and their relation to the 3D veloc- 
ity power spectrum. Section 4 shows the relation between 
observed velocity dispersion and the structure function of 
the velocity field. A way to constrain length scales (size of 
the energy containing eddies) of motions using observed pro- 
jected velocity is presented in Section 5. Method to recover 
3D velocity PS from 2D projected velocity field is discussed 
in Section 6. Discussions and conclusions are in Sections 7 
and 8 respectively. 



2 BASIC ASSUMPTIONS AND MODELS 

We consider a spherically symmetric galaxy cluster, which 
has a peaked X-ray emissivity profile. The electron number 
density is described by the /3-model profile 



ne{r) = 



no 






(1) 



where no is the electron number density in the cluster cen- 
ter (normalization) and re is the core radius. The /3-model 
provides a reasonably good description of obs erved surface 
brightness (|Cavaliere fc Fusco-Femiand 119781 ) in the cen- 
tral regions of galaxy clusters. At large radii /3— model 
is not a good descr iption of surface brightness (see e.g. 
IVikhlinin et al.ll2006l ). However the simplicity of the model 
allows us to illustrate the method and make analytical cal- 
culations. 

We have chosen /3 = 0.6 and Tc ~ 10 kpc for demonstra- 
tion of our analysis. Paramet er /3 usually varies 0.4 < /3 < 1 
for galaxy clusters (see, e.g.. IChen et al.l 120071 1. herewith a 
large fraction of clusters have /3 ~ 0.6. r^ can vary from few 
kpc to few hundreds kpc. In order to better illustrate the 
main idea of the method, we considered cool-core clusters 
with small core radius Tc (it is necessary to have a gradient 
of surface brightness down to the smallest possible projected 
distances, see Section 4 for details). 

We describe the line-of-sight component of the 3D ve- 
locity field as a Gaussian isotropic and homogeneous ran- 
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Figure 3. Expected dependence of the line-of-sight velocity dis- 
persion (line width) on the projected distance R from the cluster 
center. The increase of the effective length along the line of sight 
with the projected distance -R implies that larger and larger ed- 
dies contribute to the observed line broadening. We assume a /3 
model of the density distribution in an isothermal galaxy cluster 
(cq. [TJ and a cored 3D power law model of the velocity power 
spectrum with slope a and break at fc = km (eq. [2]l. Parameters 
are shown on the top left corner of the plot. Black dots with error 
bars: velocity dispersion with uncertainties in a single realization, 
obtained by averaging over 100 realizations (the left part of eq. 
[Sjl, red curve: analytic expression of the velocity dispersion (the 
right hand side of cq. [5] sec Appendix B). 



dom field. This allows us to gauge whether useful statistics 
could in principle be obtained. However, there is no guaran- 
tee that this assumption applies to the velocity field in real 
clusters. E.g. lEsquivel et al.l ((20071) using numerical simula- 
tions have shown that in the case of supersonic turbulence 
in the ISM, the non-Gaussianity causes some of the statisti- 
cal approaches (based on the assumption of Gaussianity) to 
fail. The same authors demonstrated that for the subsonic 
turbulence the Gaussianity assumption holds much better. 
This is encouraging since in clusters we expect mostly sub- 
sonic turbulence. Nevertheless the methods discussed here 
require numerical testing using galaxy clusters from cosmo- 
logical simulations. We defer these tests for future work. 

Power spectrum (PS) of velocity field is described by a 
cored power law model |j 



-P3D(fcx,fcy,fcz) 



B 



'^X ~l ''-"U ~i~ '^Z 



a/2 ' 



(2) 



^ Here and below we adopt the relation between a wavenumber 
k and a spatial scale x a,s k = 1/x (without a factor 27r). 
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Figure 4. Integrand from eq.|8] multiplied by k (oscillating func- 
tion) and integrand from eq. [9] multiplied by 2k (non-oscillating 
function) calculated for projected distance R = 1 kpc (black 
curves) and -R = 800 kpc from the center (red curves). The scales 

I ex — ] where these functions have maxima provide dominant 



contribution to the linc-of-sight velocity dispersion and structure 
function. The shift of the dominant scale to smaller k at large 
projected distances shows that different parts of the power spec- 
tra are probed with cr^{R) measured at different R. The shift of 
the dominant scale in the SF{l^g{R)) traces the similar shift for 
a^(R). This suggests that cr'^{R) measurements can be used as a 
proxy for the structure function at a range of scales. 



where fem is a break wavenumber (in our simple model fcm 
characterizes the injection scale) , a is a slope of the PS at 
k > km (inertial range) and B is the PS normalization, which 
is defined so that the characteristic amplitude A of velocity 
fluctuations at fc = kj-d is fixed, i.e. 

The cored power law model of the PS is a convenient de- 
scription of the PS for analytical calculations and at the 
same time resembles widely used broken power law model. 

Now let us specify the choice of parameters a and km 
in the model of the velocity PS. Cluster mergers, motions 
of galaxies and AGN feedback lead to turbulent motions 
with eddy sizes ranging from Mpc near the virial radius 
down to few tens of kpc near t he cluster core (see, e.g., 
ISunvaev. Norman, fc Brvanll2003l ). For our analysis we vary 
injection scales from 20 kpc (km — 0.05 kpc^^) to 2000 kpc 
{km = 0.0005 kpc"^). 

Parameter a - the slope of the PS - can be selected 
using standard arguments. If most of the kinetic energy is on 
large scales (injection scales), i.e. the characteristic velocity 
V decreases with k then the power spectrum PS = V^ /k^ oc 
k" with a ^ —3. At the same time the turnover time of large 



eddies should not be larger than the turnover time of small 
eddies, i.e t = l/V = l/{kV) is a decreasing function of 
the wavenumber k. Therefore, V oc k'' with 7^—1 and 
PS oc V'^ /k^ oc fc" with a ^ —5. So we expect the slope 
of the PS to be in the range —5 ^ a ^ —3. We will use 
Q = —11/3 (slope of the Kolmogorov PS), a = —4 and 
a = -4.5. 



3 OBSERVABLES AND 3D VELOCITY 
POWER SPECTRUM 

Gas motions in relaxed galaxy clusters are predom- 
inantly subsonic, and to the first approximation the 
width and centroid shift of lines measured with X-ray 
observatories contain most e ssential information on the 
ICM velocity field (see, e.g., Inogamov fc SunvaevI l2003l : 



ISunvaev. Norman, fc Brvanll2003h . That is, we have infor- 
mation about: 

(i) surface brightness in lines (Fig. (2] left panel), i.e. 
I{x,y) oc J n^{x,y,z)dz if one assumes isothermal galaxy 
cluster (effects of non-constant temperature and abundance 
of elements are discussed in Section 6) , 

(ii) emissivity-weighted projected velocity (Fig. [2] mid- 
JV{x,y,z)nl{z)dz 



die panel) V2D{x,y) 



(in practice mea- 



J n'^{z)dz 

sured projected velocity is averaged over some finite solid 
angle, see Section 5), 

(iii) emissivity-weighted velocity dispersion (Fig. [2j 
ri ght panel) a{x,y) = ^ {V'^{x,y,z));, - {V{x,y,z)y^ = 
^{V^{x,y,z)).-V^{x,y), 

where rie denotes the number electron density and V' is a 
velocity component along the line of sight. Here ()z denotes 
emissivity-weighted averaging along the line of sight, which 
we assume to be along z direction. 

Relation between the PS of the 3D velocity field, the 
2D projected velocity and the velocity dispersion for a line 
of sight with projected coordinates {x,y) are the following: 

{\V2D{kx,ky,x,y)f) = / P3o{kx,ky,k^)PEM{k:,,x,y)dk2{4:) 

and 

{cr^(a;,?y)) = / P3n{k^,ky,k^){l-PEM{kz,x,y))dkxdkydk^,(5) 



where () is the ensemble averaging over a number of real- 
izations, {\V2D{kx,ky,x,y)\^) is an expectation value of the 
2D PS of the observed projected velocity field V{x,y), Psd 
is the PS of the 3D velocity field and Pem is the PS of nor- 
malized emissivity distribution along the line of sight. For 
more details see Appendixes A and B. 

In Fig. |3] we illustrate the eq. [5] for a simple spherically 
symmetric /?-model of galaxy cluster with core radius re — 
10 kpc and /3 — 0.6. 3D velocity PS has a cored power law 
model (eq. [2}, i.e. is flat on fc < km ~ 0.005 kpc~^ and 
has a Kolmogorov slope a = —11/3 on fc > km- Observed 
velocity dispersion averaged over 100 realizations is shown 
with dots, the error bars show the expected uncertainty in 
one measurement. The right hand side of eq. [5] is shown in 
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Figure 5. Dependence of the effective length along the line of 
sight lef{{R) (eq. I12II on the projected distance R for different 
/3— models of galaxy clusters. 



red. Minor difference between two curves at small R is due 
to finite resolution of simulations. 

Once we construct the map of projected velocity 
V2u{x,y), one can also find RMS velocity of the 2D field 
at each distance r from the cluster center as 



Vnmsir) = ^{V2u{x,y)^)r - {V2B{x,y)r)\ 



(6) 



where (V2D(2;,y))r denotes mean velocity in ring at distance 
r from the center. Below we use the observed velocity dis- 
persion and RMS of the projected velocity field to constrain 
the power spectrum. 



4 STRUCTURE FUNCTION AND OBSERVED 
VELOCITY DISPERSION 

Often a structure function of the velocity field is used instead 
of the power spectrum, which is defined as 



SF{Ax) = {{V{x + Ax) - V{x)y), 



(7) 



where averaging is over a number of pairs of points in space 
separated by distance Ax. The line-of-sight velocity disper- 
sion can be linked to the structure function. Indeed, since 
the emissivity peaks at the center of the cluster and declines 
with the radius, the largest contribution to the total flux and 
to the line-of-sight velocity dispersion at distance R from the 
center comes from the region, the size of which is oc R. 

The structure function and the observed velocity dis- 
persion can be related to the PS (see appendixes C and D): 



SF{x) = 2 



and 



PiD(fcz)(l — cos 27rfcza;)dfcz 



(8) 






PlD(fcz)(l-PEM(fcz))dA;z 



(9) 



where Pid is an expectation value of the ID velocity PS and 
Pem is a PS of normalized emissivity along the line of sight. 
Fig. |4] shows the integrands in eq.[8]and eq. [9] multiplied by 
k and 2k respectively for the line of sight near the cluster 
center (black curves) and at a projected distance R = 800 
kpc from the center (red curves). Extra factor of 2 for o^{R) 
is introduced to compensate for the factor of 2 in front of 
the expression [8] for the structure function. It is more clear if 
one considers the limits of these equations at large x and R. 
When a; — >■ cxj then cos 27rfcza; oscillates with high frequency 
over relevant interval of k and mean value of 1 — cos 27rfcza; is 
~ 1. When _R — >■ cxD the emissivity distribution is very broad 
and Pem is almost a 5-function. Therefore, 



lim (SP(x)) = 2 / PiD(fcz)dfc 



and 



/oo 
Pi, 
-oo 



( yfCzjdkz 



(10) 



(11) 



From Fig. |4] it is clear that the integrands in eq. |8] and eq. 
[5] are very similar, suggesting that observed 2a^(R) should 
correlate well with the structure function. 

The structure function and the velocity dispersion (eq. 
IC5I and IDll respectively) are plotted in Fig. [S] in left col- 
umn. We fixed parameters of the /3— model of the cluster 
and varied the slope a and break km of the power spectrum 
model (eq. [2| . The relation of SF and 2cr^ is shown in the 
left bottom panel in Fig. [6] For a given R, 2a^{R) is used 
for a;— axis, while the SF is plotted as a function of lcs{R), 
where l^s is an effective length along the line of sight, which 
provides dominant contribution to the line flux, l^s is found 
from the condition that 



J^ nl{^/R?TP)dl 



:0.5. 



(12) 



Relation between ?off and projected distance depends on the 
/3— model of galaxy cluster as shown in Fig. [5] 

We then made multiple statistical realizations of the 
PS for a simple /3— model of galaxy cluster with /3 = 0.6 
and Tc — 10 kpc to estimate the uncertainties. The size 
of the box is 1 Mpc^ and resolution is 2 kpc. We assume 
that the 3D PS of the velocity field has a cored power law 
model (eq. IClf) with slope a and break wavenumber (injec- 
tion scale) at km- We made 100 realizations of a Gaussian 
field with random phases and Gaussian-distributed ampli- 
tudes in Fourier space. Taking inverse Fourier transform, we 
calculated one component of the 3D velocity field (compo- 
nent along the line of sight) in the cluster. Structure function 
and the line-of-sight velocity dispersion are evaluated using 
resulting velocity field. Right column in Fig. [6] shows veloc- 
ity dispersion along the line of sight and structure function 
averaged over 100 realizations. The expected uncertainty in 
single measurement of the velocity dispersion is shown with 
dotted curves. One can see that the overall shape and nor- 
malization of SF and 2a^ are the same as predicted from 
analytical expressions (left column in Fig.|4]), however, there 
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Figure 6. Left column: structure function (sec analytical expression IC5|I . linc-of-sight velocity dispersion multiplied by factor 2 (eq. 
IDll and eg. ID5I I and their relation for different models of cored power law power spectrum (eq. IClll . Parameters of the power spectrum 
(slope and break wavenumber) are shown in the top left corner. Parameters of the /3— model are given in the bottom right corner. 
Right column: the same as left column, but calculated by averaging over 100 statistical realizations of velocity field. The uncertainties 
in single measurement of the velocity dispersion are shown with dotted curves. 
The choice of parameters is discussed in Section 2. 
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Figure 7. Velocity dispersion and Vrms ^^ functions of distance from the cluster center with uncertainties due to stochastic nature of 
the velocity field (left panels) and their ratio from eg. 1151 fright panel). VrmsC^) ^^'^ '^{^) ^^^ averaged over a ring at distance R from 
the cluster center. Errors show the characteristic uncertainties in one measurement. The slope of assumed velocity PS is -11/3. Break 
wavenumber km of the PS is shown on the right panel. Color coding is the same in both panels. 



are minor differences (especially at small 7?) due to limited 
resolution of simulations. Relation between 2a^ and SF is in 
a good agreement with expectation relation, however the un- 
certainty in measured velocity dispersion (due to stochastic 
nature of the velocity field) is significant (see Section 6). 



5 LENGTH SCALES OF MOTIONS AND 
OBSERVED RMS VELOCITY OF 
PROJECTED VELOCITY FIELD 

Let us now consider RMS of the projected velocity field. 
During observations one gets spectra from a region, mini- 
mum size of which is set by angular resolution and/or the 
sensitivity of the instrument. RMS velocity at certain po- 
sition (x,y) for random realizations of the velocity field is 
defined as (see Appendix E) 



Cl^RMs) — 



P3D(fc.,fcy,fc,)PEM(fcz)(l-PsH(fc.,fcj;))d^fc,(13) 



where Psnikx 

fined as zero outside and unity inside the region, from which 
the spectrum is extracted (see Appendix E for details). Ve- 
locity dispersion (i.e. the line broadening) measured from 
the same region is 

(a')^ f P3u{k^,ky,k,)(l- PEMik,)PsH{k^,ky))d'k. (14) 

Therefore the ratio {VifMs)/{^^) is 

(Vims) _ I PiD{k.,ky,k,)PEM{k:,){l- Psii{k^,ky))d^k 



which can be used as an additional proxy of the length scales 
of gas motions. This ratio is mostly sensitive to the break 
of the cored power law model of the PS km- Basically at 
a given line of sight, which is characterized by an effective 
length leff the small scale motions (i.e. k > 1/Ze//) are 
mostly contributing to line broadening, while larger scale 
motions predominantly contribute to the RMS of the pro- 
jected velocity field. 

Fig. [7]shows the ratio (V^j^^g) / (a^) and its uncertainty 
calculated for different values of parameter km ■ Vrms and a 
are averaged over rings at distance R from the cluster center. 
Clearly, if km is large then all motions are on small scales 
and the ratio is small. The larger is the km, the more power 
is at large scales and the larger is the I/rms- The larger is 
the injection scale, the stronger is the increase of the ratio 
with distance and the more prominent is the peak (see Fig. 
(TJ. Here we assume that the full map of projected velocity 
field is available. Clearly, the uncertainties will increase if 
the data are available for several lines of sight, rather than 
for the full map. 

Looking at eq. [15] it is easy to predict behavior of the 
ratio on small and large projected distances R. Let us specify 
the shape of area, namely assume that we measure velocity 
in circles with radius R around the cluster center. When 
-R — >■ then the region is very small and PsH{kx,ky) — ^ 1 
over broad range of wavenumbers k < 1/R and the ratio 
Vims/f^^ -> 0. At large R Psu{kx,ky) -i> and PEM(fc^) -)■ 
(since n^iz) distribution along the line of sight on large R 
is broad), eq. 1 151 becomes 



{^^) J Psoik., ky, fc,)(l - PEM(fcz)PSH(fc., ky))d3k 



(15) 



{viMs)iR) 

{a^){R) 



J Papjkx, ky, fcz)PEM(fcz)rf^fc 
JPsod^k 



(16) 
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Figure 8. Ratio of 3D and 2D amplitudes as a function of k 
(cq. I17p for a /3-model cluster with /3 = 0.6 and re = 245 kpc 
(Coma-like cluster). A cored power law model is used for the 3D 
power spectrum, according to eq.l(2}- Two projected distances R 
were used: R = and 300 kpc. Solid lines with different colors 
correspond to different break wavenumbers in the 3D power spec- 
trum: km = 3 ■ 10^*, 3 ■ 10^'^, 3 ■ 10^^ kpc^^ and different slopes 
of the spectrum a = 3.5, 4.0, 4.5. For a given projected distance 
the curves, corresponding to different power spectra model, lay 
on top of each other, independently on the slope of the power 
spectrum, except for the case with very low break wavenumber 
km = 3 -10"^ (blue curves). For comparison the thick dashed lines 
show the ratio, calculated using simplified equation I I19I I. Clearly 
in the range of scales from Mpc and below the simple relation 
between 3D and 2D amplitudes is fully sufficient to convert the 
observed 2D power spectrum into the power spectrum of the 3D 
velocity field, unless the break of the power spectrum is at very 
large scales of few Mpc. 



and it is a decreasing function of R. 

The sensitivity of ratio Vb.ms/'^^ to the slope of the 
power spectrum is modest. There are only changes in nor- 
malization, but, the overall shape is the same. 



6 RECOVERING 3D VELOCITY POWER 
SPECTRUM FROM 2D PROJECTED 
VELOCITY FIELD 

Mapping of the projected velocity field V2D{x,y) provides 
the most direct way of estimating the 3D velocity field PS. 
The 2D and 3D PS are related according to eq|4l which we 
re-write as 



P2D{k)= f P3D{^/k^TM)PEM{k,,x,y)dk,, 
© 2009 RAS, MNRAS OOO.IlHTSl 



(17) 



where k = \/k^ + k^. This equation can be written as 
P2D{k)= I PsoiVk^ + k^,)PEM{k,,x,y)dk,+ 



oo 

+ / PsD{Vk^Tkl)PEM{k,,x,y)dk, . (18) 

l/'oft 

Contribution of the second term to the integral is small since 
PEM{kz,x,y) — >■ on fc 3> 1/leS- In the limit of fc ^ ^/^eff 
(at a given projected distance) the expression reduces to 



P2D[k) ^ PzD{k) / PEM{k„x,y)dk 



(19) 



i.e. 2D PS is essentially equal to the 3D PS of the 
velocity field apart from the normalization constant 
J PEM{kz,x,y)dkz, which is easily measured for a cluster. 
We show below that this simple relation (jlQp provides an 
excellent approximation for full expression (|17|l for a Coma- 
like clusters with flat surface brightness core. For peaked 
clusters (cool core) P30/P2U depends on projected distance 
from the cluster center since laH changes significantly with 
distance. 

It is convenient to use characteristic scale-dependent 
amplitudes of the velocity field variations, rather than the 
PS. The amplitude for 3D and 2D spectra are defined as 



A3D{k) = ^PaD{k)4TTk3 
A2D{k) = y^P2D{k)2-Kk2 



(20) 
(21) 



In these notations the relation [19] between PS transforms to 



A2D{k) = Asoik) 



1 J PEM{kz,x,y)dkz 



(22) 



Integral J PEM{kz,x,y)dkz can be estimated as - — , since 

«off 

the largest contribution is on A: < l/lctt- Eg. 1221 becomes 



A2D{k) = AsD{k) 



1 1 

2 leak 



Aaoik) 



2iVo< 



(23) 



The essence of this relation is that the amplitude of the 
3D velocity fiuctuations is attenuated in the 2D projected 



velocity field by a factor of order 



1 



A^„, 



where N^dd is the 



number of independent eddies which fit into effective length 
along the line of sight. 

We illustrate the above relation for the case of the Coma 
cluster. The density distribution in the Coma can be char- 
acterized by a /3-model with /3 — 0.6 and co-radius re — 245 
kpc. In Fig. [8] we plot the ratio A3D{k)/A2D{k) evaluated 
using equations (|19|l and (|17|l for a number of 3D PS models 
calculated at two projected distances from the Coma center. 
One can see that on spatial scales of less than 1 Mpc the 
equation (|22p is fully sufficient. The variations of the rela- 
tion for different projection distances (projected distances 
from to 300 kpc were used) affect only the normalization 
of the relation and can easily be accounted for. 

With ASTRO-H the 2D velocity field in the Coma can 
be mapped with the ~ 1.7' resolution, which corresponds to 
~ 46 kpc. Mapping 450 x 450 kpc central region of the Coma 
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would require about 36 pointings. For practical reasons it 
may be more feasible to make a sparse map (e.g. two per- 
pendicular stripes) to evaluated A2D ( e.g. cor aputing corre- 
lation function or using a method described in lArevalo et al.l 

[mi 



7 DISCUSSION 

7.1 Limiting cases of small and large scale 
motions 

Measuring characteristic amplitude of mean velocity (Vrms) 
and velocity dispersion (a) we can distinguish whether the 
turbulence is dominated by small or large scale motions. 
Clearly, the motions on scales much smaller than the effec- 
tive length along the line of sight near the cluster center can 
only contribute to the line broadening. This sets the charac- 
teristic value of the lowest spatial scale which can be mea- 
sured. The largest measurable scale is set by the maximum 
distance -Rmax from the cluster center where the line pa- 
rameters can be accurately measured without prohibitively 
long exposure time. Thus the range of scales I, which can be 
probed with these measurements is Icb{R — 0) < I < -Rmax. 
The crucial issue in measurements is "sample variance" of 
measured quantities caused by stochastic nature of turbu- 
lence. We can expect two limiting cases (see Fig. [T]). 
A: Small scale motions 

In case of small scale motions I <^ lca(R = 0) (i.e. k,n 3> 
Ich{R = 0) one expects a-{R) to be independent of radius 
and Vh,ms(^) ^ cr(_R). cr{R) is expected to have low sample 
variance and can be measured accurately even for a sin- 
gle line of sight, provided sufficient exposure time. Measure- 
ments of VRMs(fi) are strongly affected by sample variance 
and depend on the geometry of the measured map of the 
projected velocity dispersion V{x, y). If V{x, y) is measured 
only at two positions, then the uncertainty in Vrms(^) is 
of the order of its value and the ratio o{R) /Vb.ms{R) gives 
only low limit on km- We note that in this limit of small 
scale motions the assumption of a uniform and homogeneous 
Gaussian field can be relaxed and measured values of line 
broadening a simply reflect the total variance of the velocity 
along the given line of sight, while the possibility of deter- 
mining the spatial scales of motions are limited. Variations 
of a with radius will simply reflect the change of the char- 
acteristic velocity amplitude. 
B: Large scale m,otions 

In case if most of turbulent energy is associated with large 
scales (i.e. km ^ lca{R ~ 0), o-{R) is expected to increase 
with R and Vr,ms(^) ~ cr{R). In this case sample variance 
affects both (j{R) and Vr,ms(^)- Mapping the whole area 
(as opposed to measurements at few positions) would help 
to reduce the sample variance. Knowing the shape of cr(i?) 
and estimates of km from (j{R)/Vrms{R) we can constrain 
the slope of power spectrum. 



7.2 Effect of thermal broadening 

Measuring velocity dispersion one should account for line 
broadening due to thermal motions of ions. Thermal broad- 
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Figure 9. Width of the line of He-like iron line at 6.7 keV in 
Perseus cluster, assuming different ID RMS velocities. FWHM 
is related with the width as FWHM = 1.66A_B. Solid curves 
show broadening due to gas motions only, dashed curve shows 
the thermal broadening. Here we assumed the cored power law 
3D power spectrum with a = —11/3 and km = 0.001 kpc^^. 



ening should be subtracted from measured width of line. 
The broadening of the line /S.E = -\/2cr, where a is measured 
Gaussian width of the line, is defined as 



A-B 



T^ 



+ < 



(24) 



where aturb is the width due to turbulent motions and 



0"the 



fcT 

Arrir, 



is the thermal broadening for ions with 



atomic weight A. The FHWM of the corresponding line is 
FWHM = 2^/h^2A-E ^ 1.66A£;. Observed emission fines 
in the X-ray spectra of galaxy clusters correspond to heavy 
elements such as Fe, Ca, S. Because of large value of A the 
contribution of thermal broadening is small even for mod- 
est amplitudes of turbulent velocities. Indeed, the atomic 
weight for iron is 56 and the thermal width of the iron fine 
is A-B ~ 3 eV (FWHM ~ 5 eV) if one assumes the typi- 
cal temperature of clusters ~ 5 keV. At the same the gas 
motion with the sound speed would causes the shift of the 
line energy by ~40 eV. As an example, we calculated the 
expected broadening of the He-like iron line at 6.7 keV for 
the Perseus cluster, assuming Kolmogorov-like PS of the ve- 
locity field with km = 0.001 kpc~^ and varying total RMS of 
the velocity field in one dimension. The model of the Perseus 
cluster was taken from lChurazov et al.l (|2004f ) and modified 
at large distances according to Suzaku ob servations at the 
edge of the cluster (|Simionescu et al.ll201ll ). i.e. the electron 
number density is 
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Figure 10. Velocity dispersion in the Perseus cluster as function 
of distance from the center for different emissivity and abundance 
profiles. Parameters of assumed velocity PS are shown in the top 
left corner. Green: simple approximation of emissivity Ri n^; blue: 
emissivity at 6.7 keV line and flat abundance; red: emissivity at 
6.7 keV line and peaked abundance abundance (cq. I27I I in the 
center (see text for details). Black dashed curve shows expected 
uncertainty in measurements of velocity dispersion observed in a 
ring. 



7.3 The effect of radial variations of the gas 
temperature and metallicity 

The analysis described in the previous sections was done as- 
suming an isothermal /3-model spherically symmetric cluster 
with emissivity simply oc n^. Clearly, real clusters are more 
complicated, even if we keep the assumption of spherical 
symmetry. First of all the gas temperature and metallicity 
often vary with radius. These variations will be reflected in 
the weighting function n^(\/ R^ + 2^), which relates 3D ve- 
locity field and observables. To verify how strongly these 
assumptions affect the results, we calculated velocity dis- 
persion assuming detailed model of the Perseus cluster 
(see above). We assumed constant abundance profile and 
more realistic peaked ab undance profile, taken from Suzaku 
(|Simionescu et aLlbOllI ) and Chandra/XMM observations 



Z{r) = 0.4 



2-2+ f-V 



1 + 



Vso/ 



(27) 



We calculated the emissivity of the He-like iron line 
at 6.7 keV and used this emissivity as a weighting func- 
tion for the calculation of the expected projected velocity 
and velocity dispersion. Fig. [10] shows the velocity disper- 
sion along one line of sight calculated for the most simple 
model of emissivity and assuming more complicated models 
described above. One can see that mean value and uncer- 
tainties are very similar in all cases. Clearly that averaging 
velocity dispersion over a ring will decrease the uncertainty 
in one measurement (Fig. IIOI black dashed curve). 



iVe(r) 



4.68 • 10" 



1 + 



56 



+ 



4.86 • 10"^ 



i+'m)^ 



(25) 



and the temperature profile is 



1+ -^ 



T{r) = 7- 



2.3 



1 + 



5000 



(26) 



Abundance of he avy elements is assumed to be constant 0.5 
relative to Solar (| Anders fc Grevessdll9 89i). Fig.Oshows the 
calculated width of 6.7 keV line assuming various ID RMS 
velocities. Thermal broadening is shown with the dashed 
magenta curve. One can see that thermal broadening starts 
to dominate broadening due to motions only if Vrms < 150 
km/s. However, the lack of resonant scattering signatures 
in the spectrum of the Perseus cluster suggests that the ex- 
pected velocity is higher than 400 km /s in the center of the 
Perseus cluster (jChurazov et al.ll2004l ). Astro — H will have 
energy resolution of 7 eV at 6.7 keV, therefore broadening 
due to gas motions with Vrms <x 400 km/s will be easy to 
observe (Fig. [9]) in the Perseus cluster. 



7.4 Influence of density fluctuations 

The hot gas density in galaxy clusters is strongly peaked 
towards the centre. However, besides the main trend, there 
are density fluctuations, which could contribute to the ob- 
served fluctuations of the projected velocity. Let us split the 
density field in two components 



n — no + Sn, 



(28) 



where no corresponds to the smooth global profile and Sn 
represents fluctuating part. Neglecting terms of the order 
(5n) , the emissivity-weighted projected velocity is 



V2ii{x,y) 



J no{z)V{x, y, z)dz + 2 J no{z)SnV{x,y, z)dz 
Jnl{z)dz 



(29) 



The analysis of X-r ay surface brightness fluctuations in the 
Coma cluster (e.g. IChurazov et al.ll2012l 'l have shown that 
the amplitude of the density fluctuations ( — ) is of the order 
~ 2 — 10%. Assuming that the same number is applicable 
to other clusters, the second term in eq.[29]is small and the 
contribution of the density fluctuation to projected velocity 
can be neglected. 



7.5 Measurements requirements 

To illustrate the most basic requirements for the instru- 
ments to measure the ICM velocity field, let us consider 
two examples of a rich cluster and an individual elliptical 
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Table 1. Energy and angular resolutions requirements for future 
X-ray observatories. 

Perseus eluster NGC5813 
FeXXV line OVIII line 



FWHM (therm.) 

Shift of line energy (V= 

Angular resolution 



:e.) 



4.9 eV 


0.32 eV 


25.7 eV 


0.91 eV 


0.5' 


2" 



galaxy, representing the low temperature end of the cluster- 
group-galaxy sequence: 



(i) Perseus cluster, Tmean 
keV, distance 72 Mpc 

(ii) Elliptical galaxy NGC5813, T^ea 
OVIII at 0.654 keV, distance 32.2 Mpc 



= 5 keV, line of FeXXV at 6.7 



0.65 keV, line of 



Table [T] shows the FWHM of these lines calculated if 
only thermal broadening is taken into account, shift of the 
line energy in the case of gas motions with V — Cs and 
desirable angular resolution of the instrument. Future X- 
ray missions, such as Astro-H and ATHENA will have en- 
ergy resolution ~ 5 eV and ~ 3 eV respectively. Such en- 
ergy resolution is sufficient to measure broadening of lines in 
hot systems like galaxy clusters. Cold systems, like ellipti- 
cal galaxies, require even better energy resolution. However, 
turbulence still can be measures using the resonant scat- 
terin g effect (see e.g. IWerner et al.ll2009l : IZhuravleva et al.l 
I2OIII ) or by using grating spectrometers observations. An- 
gular resolution of ~ 1' for nearby galaxy clusters should 
sufficient to study the most basic characteristics of the ICM 
velocity field, while for elliptical galaxies the resolution at 
the level of ~ arcsec (comparable with Chandra resolution) 
would be needed. 

Astro-H observatory will have energy resolution ~ 5 eV, 
field-of-view 2.85' and angular resolution 1.7', which means 
that it will be possible to measure shift and broadening of 
lines as a function of projected distance from the center 
in nearby clusters. E.g. if one assumes that RMS velocity 
of gas motions in Perseus cluster is ~ 300 km/s, then ~ 
4 ■ 10^ s is enough to measure profiles of mean velocity and 
velocity dispersion with a statistical uncertainty of ~ 30 
km/s (90% confidence) in a stripe ± 200 kpc (7 independent 
pointing and 28 independent measurements in 1.4' x 1.4' 
pixels) centered in the center of the cluster q In order to 
measure velocity with the same accuracy at larger distances 
from the center , e.g. at 500 kpc and 1 Mpc, one would need 
8 • 10^ s and 4 • lO*" s exposure respectively. 



^ The estimates were done using the cur- 

rent version of Astro-H response at 

|http://astro-h. isas.jaxa.jp/resear chers/sim/response.html| 



8 CONCLUSIONS 

Various methods of constraining the velocity power spec- 
trum through the observed shift of line centroid and line 
broadening are discussed. 

• Changes of the line broadening with projected distance 
reflects the increase of the spread in the velocities with dis- 
tance, closely resembling the behavior of the structure func- 
tion of the velocity field. 

• Another useful quantity is the ratio of the characteristic 
amplitude of the projected velocity field to the line broad- 
ening. Since the projected velocity field mainly depends on 
large scale motions, while the line broadening is more sensi- 
tive to small scale motions, this ratio is a useful diagnostics 
of the shape of the 3D velocity field power spectrum. 

• Projected 2D velocity field power spectrum can be eas- 
ily converted into 3D power spectrum. This conversion is 
especially simple for cluster with an extended flat core in 
the surface brightness (like Coma cluster). 

Analytical expressions are derived for a /3-model clusters, 
assuming homogeneous isotropic Gaussian 3D velocity field. 
The importance of the sample variance, caused by the 
stochastic nature of the turbulence, for the observables is 
emphasized. 
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APPENDIX A: 3D VELOCITY POWER 
SPECTRUM AND PROJECTED VELOCITY 
FIELD 

Let us assume that the line-of-sight component of the 3D ve- 
locity field V3D(a;, y, z) is described by a Gaussian (isotropic 
and homogeneous) random field. We assume that the cen- 
troid shift and the width of lines contain most essential in- 
formation on the velocity field. Here and below we adopt 
the relation k — 1/r without factor 27r (see Section 2 for 
details). 

Projected 2D velocity along the line of sight (observed 
centroid shift of the emission line) in z direction is 



V2D{x,y) 



J V3d{x, y,z)nl(x,y, z)dz 
Jni{x,y,z)dz 



(Al) 



where Vsd is z component of the 3D velocity field and no is 
the electron number density. Denoting normalized emissivity 
along the line of sight at a certain position with coordinates 
(x,y) as e{z) = n^{x,y,z)/ J n^{x,y,z)dz the previous rela- 
tion can be re-written as 



V2o{x,y)= I V3n{x, y,z)e{z)dz. 



(A2) 



Applying the convolution theorem one can find the Fourier 
transform of Vsu{x,y, z)e{z) as 



V3o{kx,ky,k2-^)e{k^ - k^^)dk^-^, 



(A3) 



where Vsd and e are Fourier transforms of the 3D velocity 
field and normalized emissivity respectively. The projection- 
slice theorem states that 



f2D{kx,ky) ^ f3Y){kx,ky,0). 



(A4) 



Accounting for IA3I and IA4I we can write the Fourier trans- 
form of 2D velocity field as 



V2D{ka:,ky) = / Vsoika: , ky , kz-^)t {kzj^)dkzj^ , 



(A5) 
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where * denotes conjugation. 

Averaging over a number of realization we find the 
power spectrum of the projected mean velocity 

{|V2D(fc.,fcJ|') = (j j V3u{k^,ky,k,,)t{k,,)dk,,\'). (A6) 

The right part of the equation above can be re-written as 

{Vw(k:c,ky,k^^)V-^-o{kx,ky,k^.^))x 

xt {k^^)e{k^2)dk^^dk^2. (A7) 

Since phases are random, all cross terms after averaging over 
a number of realizations will give if fc^j 7^ kz^- Therefore, 

{|V^2D(fc.,fcj,)|')= f\V3o{k^,ky,k,,)f\e{k,,)fdk,,. (A8) 

Denoting power spectra of the 3D velocity field and nor- 
malized emissivity as P3D and Pem respectively the final 
expression is 

{\V2o{k,,ky)f) = [ PMkz,ky,k,)PEM{k,)dk,. (A9) 



APPENDIX B: 3D VELOCITY POWER 
SPECTRUM AND PROJECTED VELOCITY 
DISPERSION 

Projected mean velocity dispersion for the line of sight with 
coordinates (x,y) averaged over a number of realization is 
defined as 



(o- (s,y)) = ( / V3uix,y,z)e{z)dz)' 



-{( / V3oix,y,z)e{z)dzj ). 

It can be re-written as 

(o"^(a:,y))= {V3-a(x,y,z))e{z)dz- 



(Bl) 



-{Viu{x,y)}. 



(B2) 



Expanding V3D and V2D through the Fourier series, averag- 
ing over realizations and keeping non-zero cross terms will 
give 



/ 



/• 



{a {x,y)) = / \V3n{kx,ky,kz)\ dk^dkydk^ / e{z)dz 



\V2vikx,ky)\ dkxdky. (B3) 



Since emissivity along the line of sight is normalized so that 
J e{z)dz = 1 and accounting for ea lA9l the final expression 
for projected velocity dispersion is 

{a\x,y)) = f P3uik^,ky,k,)(l~ PEMik,))dk:,dkydk,.(B4) 



APPENDIX C: RELATION BETWEEN 
STRUCTURE FUNCTION AND CORED 
POWER LAW 3D POWER SPECTRUM 

Let assume 3D isotropic and homogeneous power spectrum 
(PS) of the velocity field is described as 

B 



P[io(k^,ky,k^) = 






(CI) 



where fcm is a wavenumber where /3 model has a break (e.g. 
an injection scale in turbulence model), a is a slope of PS 
on fc > fcm and B is PS normalization, which is defined so 
that the characteristic amplitude A of velocity fiuctuations 
at fcrof is fixed, i.e. 

B = . ., „ ,. r. (C2) 



47rfc3^fP3D(fcrcf) 



Integrating ICll over dk^dky = 2-Kkdk, one can find ID 



PS 



Pid(A:z 



2 \ -"/2+1 
2^Bkl ( 1 + ^ 

OL-2 



(C3) 



Structu re function (SF) is related to ID PS by the transfor- 
mation iJRvtov et al.lll988l l 

/-t-00 
PiD(fcz)(l-cos27rfea;)dfcz. (C4) 

-00 

Substituting IC3I to IC4I and assuming a > 3 yields 



SF{x) 



4Bki (7r^r(C) -27r*fcia;«i^c(C,27rsfcnO) 
(a-2)r(f-l) 



(C5) 



where ?=§■ — § and Ka is a modified Bessel function of 
the second kind. 



APPENDIX D: RELATION BETWEEN 
VELOCITY DISPERSION ALONG THE LINE 
OF SIGHT AND POWER SPECTRUM 

PiD and P3D are related as Pio{kz) = 

J P3u{kx,ky,kz)dkxdky, therefore eg. IB4l can be re-written 
as 



(a^R)) = / PlD(fcz)(l - PEM(fc.))dfc. 



(Dl) 



where R is projected distance from the center. 

If electron number density is described by /3 model with 
normalization no and core radius Re 

no 



na{r) 



1 + 



R^ + z' 
RI 



(D2) 



then the emissivity is 

where r^ = R? -\- z"^ , C = R^ + R? and we assume normal- 
ization 710 = 1. The Fourier transform of emissivity is 



(D3) 
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J —(. 



■ cos{2'Kkzx)dx^ 



(D4) 



(C7 + a;2)3'3 

where terms with isin{2-Kkzz) axe zero since integrand is 
symmetrical. Dividing ID4I by the total flux and assuming 
P > 1/6 one can find weight Vy(fcz) as 



W{k, 



2cHi-K^K^{-C,2VCK-R) 

r(C) 



(D5) 



where (^ = 3/3 — i and Ka is a modified Bessel function of 
the second kind. 



APPENDIX E: RATIO OF OBSERVED RMS 
VELOCITY TO OBSERVED VELOCITY 
DISPERSION 

Let us assume that spectrum is extracted from the region 
with the area J dxdy. The RMS velocity from this spec- 

shape 

trum averaged over a number of realizations is 



(^RMs) 



/ V^j^{x,y)dxdy 

'haps 

J dxdy 

J V2r>(.'^,y)dxdy\ 
f dxdy 



(El) 



where () denotes averaging over realizations. Accounting for 
ea. lA2l we can re- write the first term in the above equation 
as 

/ (SV3D{x,y,z)e(z)dzfdxdy 
shape 



J dxdy 



"Td^y I /<^3D(fc.„fc,i,fc.J)e''"'===i^e"^'="^x 



shape 



shape 



X Fem (fc^i ) { VsD {kx2 , ky2 ,kz^))t 



i27rfca:2 ^^^^Trfcy^y 



-FEM(fcz2)x 



Xtfkicfk2dxdy, 



(E2) 



where -Fem is a Fourier transform of emissivity along the line 
of sight. Averaging over a number of realizations will leave 
non-zero terms only if ki — k2. Therefore, the first term in 
eg. [Ellis 

/ J P3o{kx,ky,k,)PEM{kz)d''kdxdy 

shape 

J dxdy 

shape 

= J PaD{kx,ky,kz)PEM{kz)dkxdkydkz. (E3) 

The second term in eq. lEll can be written as 



/ f V3-D(x,y,z)'i(z)dxdydz\ 
f dxdy 



r ^ — *27rfcT.. x-i „ — z27rfc,,, V] vy 

X J e °=i ^e VI"'- X 

shape 
^ J e^'^^h^^X2^i2,,ky^y2^^_^^y^^^^^y^^a^^ 

shape 



(E5) 



/ dxdy j 
'.hapc J 



-I J P3Dikx,ky,kz)PEM{kz)x 



X / 5(2:2, y2)e'^''''-^''^e'^''''y2y^dxidyidx2dy2d^k, 



(E6) 



where S{x, y) is a mask, which is zero outside and unity 
inside the area, from which spectrum is extracted. Denoting 
power spectrum of 5(2;, j/)/( / da::dj/)^ as PsH(fci, fcy), the 

shape 

second term in eg. lEll becomes 



Pso{kx, ky, k,)PEM{kz)PsH{kx, ky)dkxdkydkz 



(E7) 



Therefore, the final expression for the RMS of the projected 
velocity field is 

(V^RMS) = f P3u{kx,ky,k,)PEM(k,)(l~PsH{kx,ky))d^k.{E8) 

And the ratio of VjIms ^nd velocity dispersion along the line 
of sight is 

(Vims) _ I Pw{kx,ky,k,)PEM{k,){l~ PsH{kx,ky))d''k 



(a^) 



S Pwikx,ky,k,){l- PEM{k,))d^k 



(E9) 



This paper has been typeset from a T^hJX/ I^TJhJX file prepared 
by the author. 



(E4) 



/ /V'3D(fcx,feB>fe.)<='^""'°"e'^"'''° 



BFemCsj:)'* *:da:dt/ 



r dxdy 



Squaring and averaging over realizations yield 

^ J P3i){kx,ky,kz)PEM{kz)x 



f dxdy 
■■hape 
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